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Abstract 

A generalized screened potential model (GSPM), recently developed to stndy the bot- 
tomoninm spectrnm, is applied to the calcnlation of charmoninm masses and electromag¬ 
netic widths. The presence in the GSPM of more qnark-antiqnark bonnd states than in 
conventional non screened potential models, allows for the assignment of GSPM states to 
cataloged non conventional charmoninm resonances as well as for the prediction of 
new (non cataloged) states. The resnlts obtained seem to indicate that a reasonable 
overall description of charmoninm resonances is feasible. 
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1 Introduction 


In a recent paper [T] a new non relativistic quark model for the description of heavy 
quark mesons has been developed. The novelty of the model, called Generalized Screened 
Potential Model or GSPM, is the consideration of a lattice motivated quark-antiquark 
interaction that implicitly incorporates color screening effects from meson-meson conhg- 
urations. When applied to bottomonium a good spectral description of well established 
resonances is obtained and a richer high energy spectrum (bigger number of bound states) 
than the one resulting from the non-screened Gornell potential is predicted. However, the 
current lack of data does not allow to validate or refute this prediction. In this regard the 
application of the model to charmonium could be determinant since a plethora of addi¬ 
tional states, not htting into the conventional non-screened Gornell potential framework, 
has been discovered in the last ten years (see PEIS] and references therein). 

In this article we apply the GSPM to charmonium. We extend the observable analysis 
beyond the spectral masses to electromagnetic widths, for the model is suitable for their 
calculation and there exist data to be compared with. We do not analyze strong decays 
since a fully consistent treatment of them within the GSPM framework (involving the 
description of mesons containing light quarks) is a formidable task outside the scope of 
the present study. 

We show that a reasonable description of well established and candidates to 
resonances is feasible. Moreover the model allows for some dehnite predictions about 
new resonances what might be used in future experimental searches to further check its 
validity. Regarding 1 states the presence of overlapping thresholds limits the appli¬ 
cability of the GSPM to spectral energies quite below the hrst meson-meson threshold 
becoming then completely equivalent to the Gornell model. 

The article contents are organized as follows. In Section [2] a brief review of the GSPM 
is presented. In Section [3] the model is applied to the calculation of the charmonium 
spectrum and electromagnetic widths and the results are compared to data. A calculation 
from a non screened Gornell potential is also shown for comparison. Finally, in Section 
[Hour main results and conclusions are summarized. 

2 Generalized Screened Potential Model (GSPM) 

The Generalized Screened Potential Model (GSPM) is based on the assumption that a 
heavy quark meson description can be attained from the consideration of effective valence 
quark degrees of freedom interacting through a potential that incorporates screening 
effects from meson-meson conhgurations. 

More precisely, the Generalized Screened Potential that we shall call V (r) henceforth 
tries to implement within a quark model framework the lattice results for the energy of 
two static color sources (heavy quark and heavy antiquark) in terms of their distance, 
Eiattice (’"), when the mixing of the quenched quark-antiquark conhguration with open 
flavor meson-meson ones is taken into account. In reference [S] the lattice calculation 
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for the case of one open flavor meson-meson configuration was performed, the resulting 
Eiattice (^) having a different form below and above the meson-meson threshold. For the 
two threshold case an educated guess for Eiatuce i.'r) was done (see Fig. 22 in [5]). A 
simplihed generalization of these lattice results to the many threshold case was proposed 
in reference [1] . From it a static quark-antiquark potential, V (r), was derived by means 
of a Born-Oppenheimer approximation, say by subtracting the quark and antiquark 
masses, tuq and m-g, from the static energy. 

Explicitly, by calling MTi with i > 1 the masses of the physical meson-meson thresh¬ 
olds, Ti, with a given set of quantum numbers and dehning Mtq = 0 for a unihed 

notation (note that Tq does not correspond to any physical meson-meson threshold), the 
form of V (r) in the different energy regions (specihed as energy interval subindices) 
reads: 
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with the parameters a and y standing for the string tension and the color coulomb 
strength respectively. 

Thus V (r) has in each energy region between neighbor thresholds a Cornell form but 
modulated at short and long distances by these thresholds. 


Thus for example in Fig. [T]the form of V (r) in the hrst and second energy regions 
is drawn for cc states with = 0+(l’''’'') quantum numbers, whose hrst threshold 

Ti corresponds to DD* and its second threshold T 2 to D*~. 

Let us remark that V (r) is an energy dependent potential in the sense that its form 
differs in the different energy regions delimited by the thresholds. Actually this is the 
essential difference with other screened potential models which have been also employed 
for the description of charmonium 

It is also important to emphasize that in any energy region the potential is strictly 
conhning in the sense that only bound states are obtained as solutions of the Schrodinger 
equation in such energy region (see next Section). 


2 





Figure 1: Generalized screened potential V{r). The solid (dashed) line indicates the 
potential in the hrst (second) energy region for 0’''(1’''’'') cc states with rric = 1348.6 
MeV, a = 850 MeV/fm, y = 100 MeV.fm, Mt^ = 3872 MeV and Mt^ = 4080 MeV 
(values of the parameters and threshold masses from Section [3l). 


3 Charmonium 

From the dehned V (r), charmonium states in the energy region , char¬ 
acterized by a dehnite set of quantum numbers are obtained by solving the 

Schrodinger equation for 

In order to get the solutions we previously £x the values of the parameters of the 
model and list the open charm meson-meson threshold masses to be considered. Then 
we detail the calculation of the spectrum for a particular case before giving the gen¬ 
eral results. Next we assign calculated states to experimental resonances and use the 
corresponding wave functions to evaluate electromagnetic widths. 

3.1 Parameters 

Let us realize that for energies quite below the first corresponding thresholds the potential 
V (r) is almost completely equivalent to a Cornell one: 

Vcor{r) = or - - (r : 0 CX)) (4) 

r 

But the use of the conventional Cornell potential teaches us that the charm-anticharm 
(cc) system may be a relativistic one [9]. This makes debatable the application of the 
GSPM to charmonium. In the spirit of quark model calculations we shall assume that the 
effectiveness of the parameters (quark mass, string tension and coulomb strength) may be 
appropriately taking into account, at least in part, relativistic corrections. We can also 
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invoke the effectiveness of the parameters regarding additional contributions from light 
quark-antiquark pairs apart from the implicitly considered meson-meson configurations. 

Let us also note that i) no threshold widths have been considered and ii) the ac¬ 
cumulative interacting effect from different meson-meson configurations with the same 
threshold mass has not been implemented. (Indeed if a meson-meson configuration gives 
rise to a screening of the quark and antiquark color charges for a given energy then a 
reinforcement of the screening is expected when more meson-meson configurations with 
the same threshold mass are available.) 

Hence we shall restrict the application of the model to those energy regions involving 
non degenerate, isolated (in the sense of not having any significant experimental overlap 
due to their widths) thresholds. 

Even so the model may be too simplistic for an accurate description of real mesons. 
On the one hand V{r) does not contain spin dependent terms that we know may give 
significant contributions to the masses of the lower spectral states (see for example [IQ]). 
On the other hand the effect of any threshold has been approximated by an abrupt 
(instead of a physically soft) change in the potential at the crossing radii. Moreover 
SU (3) flavor symmetry has been considered when the same effect (flattening of the 
potential from the crossing radii) from thresholds with ss, uu or dd content has been 
implemented despite the fact that the probability of formation for each of these pairs 
may be different. 

Keeping in mind these possible shortcomings we shall try to show that such a sim¬ 
ple model could provide us with some insight onto the dominant dynamic mechanisms 
governing the charmonium structure. 

As we are dealing with a spin independent potential we shall compare as usual the 
calculated s— wave state masses with spin-triplet data, the p— wave state masses with the 
centroids obtained from data and the d— wave states with the few existing experimental 
candidates. 

Aiming at a joint description of charmonium and bottomonium we shall use for both 
the same values for the parameters of the potential. From [1] we have a = 850 MeV/fm 
and X = 100 MeV.fm. Let us realize that this string tension value y/a = 410 MeV is 
within the interval usually accepted for it from phenomenology (see for instance [H]). As 
for the Coulomb strength y its value corresponds to a strong quark-gluon-quark coupling 
Os = ^ ~ 0.38, in agreement with the value derived from QCD from the fine structure 
splitting of Ip states in charmonium [12]. Regarding the remaining parameter of the 
model rric we fix its value to get a reasonable overall fit to the spectrum. 

Thus the set of parameters that will be used henceforth is 

a = 850 MeV/fm 

X = 100 MeV.fm (5) 

TTic = 1348.6 MeV 

where the value of the charm mass rric = 1348.6 MeV has been fine tuned to fit the mass 
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of the well established non conventional charmonium state X(3872) since this resonance 
may be naturally described in the GSPM as explained later on. 

It is important to remark that due to the effective character of the parameters a 
better overall spectral £t (differences from calculated masses to data of 35 MeV at most) 
could be achieved by choosing for example a = 750 MeV/fm, y = 100 MeV.fm and 
rric = 1407.8 MeV. However as this new fit makes no difference at all in the resulting 
number of spectral states we prefer to maintain the same potential description as in 
bottomonium. 

3.2 Thresholds 

In order to apply the GSPM to a particular set of charmonium (/ = 0) states with 
definite we need the masses for open charm meson - meson thresholds coupling 
to these quantum numbers. From these masses the crossing radii are immediately 
calculated from ([3]). 

Let us realize that the static approach we follow to build the potential implies that 
the heavy quark and antiquark or, quite equivalently, the two charmed mesons forming 
the threshold, are in a relative S— wave so that the threshold mass is just the sum of 
the masses of the mesons. 

The list of known thresholds, their masses and the corresponding crossing radii 
appear in Tables [T] and [2] where a simplified notation has been used: a threshold 
has been denoted by the hrst meson-meson component entering in the / {J^^) lin¬ 
ear combination. Thus, the first threshold in Table [H D*^(2007) denotes 

(L)OI^(2007), D+D*(2010)-) 1^0; jp=i+ +C.C. where c.c. stands for the charge conjugate. 


We have used isospin symmetry to construct thresholds with well dehned isospin. 
This means that we are neglecting the mass differences between the electrically neutral 
and charged members of the same isospin multiplet, for example and with PDG 
quoted masses [2] 1864.91T0.17 and 1869.5±0.4 respectively or i7*(2007)° and Zl*(2010)“ 
with quoted masses 2006.98 ± 0.15 and 2010.21 ± 0.13 respectively. For the calculation 
of the threshold masses we have used the lower mass value in any isospin multiplet (1865 
MeV and 2007 MeV in the examples just mentioned). 

Regarding the C parity, for a threshold formed by two mesons 91li and VJI 2 we can 
construct the combinations (9Jti9Jt2 i c.c.) with C parity -|- and — respectively. Notice 
though that if 2 JI 2 = then, as the two mesons are in S— wave, we have DJliDJli = 
(^_yi+n-3 where ji stands for the spin of and j for the total spin of the 

threshold. Therefore only one combination in c.c. is allowed for a given value 

of j (the other vanishes). For example the J = 0 threshold D*D* with ji = 1 has only 
positive C parity when coupled to j = 0, 2 (and only negative C parity when coupled to 
J = !)• 
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Ti Charmonium Thresholds (MeV) (fm 

0(0++) 


Ti 


(o-,o-) 

3730 

1.31 

T 2 

DtD- 

( 0 -. 0 -) 

3937 

1.54 

Ts 

T)*°(2007)I^(2007) 

( 1 -. 1 -) 

4014 

1.62 

T 4 

D*+D*- 

(1-.1-) 

4224 

1.86 

T 5 

(2550) 

(o-,o-) 

4405 

2.07 


0(1++) 


Ti 

£>”£>• (2007) 

(o-,i-) 

3872 

1.46 

T 2 

DtDT 

(0-.1-) 

4080 

1.70 


0(2++) 


Ti 

-D*(2007)‘>D*“(2007) 

(i-.iO 

4014 

1.62 

T 2 


(1-.1-) 

4224 

1.86 


Table 1: Open charm meson-meson thresholds for 0 (J++) charmonium states. and 
J 2 stand for the angular momenta of the mesons forming the threshold. Threshold 
masses {Mt^) obtained from the charmed and charmed strange meson masses quoted in 
|2]. Crossing radii (r^J calculated from ([3]). 
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I (J'’°) 

Ti 

Charmonium Thresholds 

( Jf , ) 

Mxi 

(MeV) 

Tt- 

(fm) 

0(1--) 

Ti 

D0I77°(2420) 

L)Od7°(2430) 

(o-,i+) 

4287 

1.93 


T 2 

i:)*(2007)°D*°(2400) 

(i-,o+) 

4325 

1.98 



T>*(2007)°I5^(2420) 






T>*(2007)°Dl^(2430) 


4429 

2.09 


-*■ 0 

D+D,i(2460)- 

(o-,i+) 



D:+D:o(2317)- 

(i-,o+) 





D*{2007)°Df{2AQ0) 

(l-,2+) 

4470 

2.14 


Ts 

D+D,i(2536)- 

(o-,i+) 

4504 

2.18 


T, 

D*+Dsi{2m)- 


4572 

2.26 


Tj 

D:+T),i(2536)- 


4648 

2.35 


Ts 

D:+D:2(2573)- 

(l-,2+) 

4685 

2.39 


Table 2: Open charm meson-meson thresholds for 0 (1 ) charmonium states. Ji and 

J 2 stand for the angular momenta of the mesons forming the threshold. Threshold 
masses obtained from the charmed and charmed strange meson masses quoted in 

[2]. Crossing radii (rjv) calculated from ([3]). 
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A look at Table [U makes clear that for 0 states only non degenerate isolated 

thresholds are present. Therefore the GSPM can be safely applied. On the contrary for 
0 (1 ) states, Table 121 there are degenerate (Ti and T 3 ) and overlapping thresholds (for 

instance Ti and T 3 overlap with T 2 due to the large width (267 MeV) of 15^(2400)). 
Therefore, at its present stage the GSPM can only be consistently applied to 0 (1 ) 

states quite below the hrst threshold. In consequence we shall restrict our study in this 
case to the hrst energy region. 


3.3 Spectrum 


Gharmonium (cc) states are obtained by solving the Schrodinger equation for V{r). As 
in any energy region it is a radial potential we use spectroscopic notation k = nl, in 
terms of the radial, n, and orbital angular momentum, I, quantum numbers, to denote 
its bound states. Thus in the energy region we have 


= . ^,1 (cc)fc 






( 6 ) 


where T stands for the kinetic energy operator, 
, for their masses. 



""[Ti-l.Ti] 


for the bound states and 


Let us consider for example the 0’''(1’''’'') spectral states. In the hrst energy region 
the potential V{r) has the form j(r), given by ([ID (solid line in Fig. HD 

( <jr — ^ r < 1.46 fm 

^[ 0 , 3872 ] (?") = s (7) 

[ 1174.8 MeV r > 1.46 fm 

where and have been taken from Table [1] and the values of the parameters 
(a, X, rric) are given by ([5D. 


By solving the Schrodinger equation for Vio, 3872 ](’") we get the GSPM 0’''(1+’'') spec¬ 
trum in the hrst energy region [Mt^ = 0, Mti = 3872 MeV]. It has two bound states 
states, lp]^To,Ti] and 2 p[ 7 ’g^rp, whose masses generically denoted by Mgspm are 

listed in Table O 


In the second energy region, [Mt^ = 3872 MeV, Mt^ = 4080 MeV ], the potential 
V{r) has the form j(r), given by (|2D (dashed line in Fig. [ID: 


' 1174.8 MeV r < 1.46 fm 

V[ 3872 , 4080 ](^^) = < OT - ^ 1.46 fm < r < 1.70 fm 

1382.8 MeV r > 1.70 fm 


( 8 ) 
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cc 

0(1++) 


[dTp_i, VfrJ 
MeV 

GSPM 

States 

hTi-i,Ti] 

Mgspm 

MeV 


[To,Ti] 

[0, 3872] 

lP[7o,ri] 

3456.1 




2p[ro,ri] 

3871.7 


[Ti,T,] 

[3872,4080] 

IPTi.Ta] 

4017.3 


Table 3: Calculated charmonium masses from V{r), generically denoted by 

Mgspm, in the first two energy regions specified by the thresholds [Tj_i,Tj] and their 
masses 

where the threshold masses and crossing radii are taken from Table [H The spectrum 
has only one bound state lp[Ti,T2] whose mass generically denoted by Mgspm 

is listed in Table [3l 


By proceeding in the same way for higher energy regions and for different quantum 
numbers we get the complete GSPM bound state spectrum. But before listing the 
calculated spectral masses it may be illustrative to analyze the effect produced by just 
one threshold. For this purpose we shall compare the results obtained from the Cornell 
potential with the ones obtained from a GSPM with only one threshold. 

Let us consider again states and calculate the spectrum if only the thresh¬ 

old Ti (corresponding to D^D* (2007)) is present. Then there will be two energy re¬ 
gions. In the hrst one, [Mtq = 0, Mt^ = 3872 MeV], the potential is given by ([7]). Hence 
there are two bound states, lp^To,Ti] and 2p[ro with the masses previously calculated 
(see Table [3]) which have been listed again in Table HI In the second energy region, 
[Mti = 3872MeI/, oo[, the potential reads 

( 1174.8 MeV r < 1.46 fm 

hi3872,oo[(’") = \ 

[ err — ^ r > 1.46 fm 

and has an inhnite number of bound states. The masses of the two lowest states in this 
energy region are listed in Table HI 


For the sake of comparison we calculate the Gornell spectrum in the same energy 
interval (from 0 to 4400 MeV) from the same values of the parameters a, x nnd rric 
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GSPM 

States 

Mgspm 

MeV 

Mgot 

MeV 

Cornell 

States 

lP[7o,ri] 

3456.1 

3456.2 

Ip 

2P[ro,ri] 

3871.7 

3910.9 

2p 

lp[ri,oo[ 

4029.3 



2p[Ti,oo[ 

4303.3 

4294.6 

3p 


Table 4: Calculated charmonium masses up to 4400 MeV when only the first 

threshold is considered: Mgspm- Masses from the Cornell potential, Mcor-, are also 
shown for comparison. Conventional spectroscopic notation has been used to denote the 
Cornell states. 

given by ([5]). The results are also listed in Table HI We immediately realize that in the 
spectral region considered there are four GSPM states [lp[To,Ti], ‘^Pito,Ti], lp[ri,oo[, ‘^P[Ti,oo[) 
for only three Cornell states {lp,2p,3p). Moreover, the masses of the hrst and fourth 
GSPM states are quite the same as the masses of the hrst and third Cornell states. Since 
the GSPM potential differs from the Cornell one in the incorporation of the threshold 
Ti = D^D* (2007) we may interpret these results by saying that the second and 

third (lp[Ti,oo[) GSPM states are effectively describing the mixing of the second Cornell 

(2p) state with the D*^(2007) conhguration. Therefore the effect of the threshold is 
the appearance of one more spectral state (notice though that if the 2p Cornell state 
were farther above the threshold the GSPM would not generate the 2pp-f^^Ti] state). 

It may also be interesting to compare the resulting radial wave functions for the 
2P[To,Ti] GSPM state and the 2p Cornell state. This comparison is drawn in Fig. |2l 

As can be checked the 2p[To,Ti] radial wave function extends to much larger distances 
than the 2p one. If we consider the 2p[Xo,Ti] state as an effective description of the 
experimental X(3872) and the 2p Cornell state as describing a (non experimental) con¬ 
ventional Xci (2p) state then it is clear the difference between them. The comparison of 
the respective root square mean radii, 3.6 fm for X(3872) and 1.1 fm for Xci (2p), indi¬ 
cates the screening of the heavy quark color charges in X(3872) due to the presence of 
the threshold. (Let us point out that in a couple channel treatment involving quenched 
quark-antiquark and meson-meson conhgurations this would correspond to the presence 
of a D^D* (2007) wave function component.) 
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Figure 2: Radial wave functions R(r) (in units /m“i) for the (2p^To,T{\) GSPM state 
(thick line) and the (2p) Cornell state (thin line). 



3.3.1 0+ (J++) GSPM States 

The spectrum for 0^(J’''+) cc states from V{r) is shown in Table O The spectrum from 
the Cornell potential Vcorif) given by ([1]) with the same values of the parameters a, x 
and rric is also listed for comparison. 


A glance at the table confirms the presence of a bigger number of CSPM states than 
Cornell ones even ignoring possible additional I"*"’'' states above 4080 MeV and 2~^~^ states 
above 4224 MeV. More precisely there are (at least) four J~^~^ CSPM states in the energy 
interval 4000 — 4400 MeV for only one Cornell state. Since the calculated masses of three 
of these CSPM states are in good correspondence with the masses of X (4140) O’*" 
and X (4350) O’*" (?’"''), the currently existing experimental candidates to 0'''( J'*""'') states 
in that energy interval (see Particle Listing in [2]; see also 0), a tentative assignment 
of CSPM states to these candidates has been done in Table O From it a guess for their 
unknown quantum numbers comes out: X (4140) O’*" (O’*"’*') or X (4140) O’*" (2’''’'') and 
V (4350) 0+ (0++). 

Furthermore the model predicts the existence of at least two new 0’''( J’*"*') resonances 
in the energy interval considered. One of them, that we shall call 0(4140) {C stand¬ 
ing for theoretical candidate) would be assigned to the 2~^~^ 0++ (lP[T3,r4]) 

CSPM state at 4140.2 MeV (see Tabled]). Let us note that the existence of this state 
in the CSPM is linked to the existence of the 0’''+ (lP[r 3 ,r 4 ]) or 2++ (1 P[Ti,T 2]) state 
that we have assigned to X (4140); as both states are in between the same thresholds 


11 




JPC 


GSPM 

States 

hTi-i,Ti\ 


Mqspm 

MeV 


MpDG 

MeV 


Mcor 

MeV 


Cornell 

States 

k 


0++ 

1++ 

2++ 

lP[To,Ti] 

IPTo.Ti] 

IPTo.Ti] 

3456.1 

3456.1 

3456.1 

3414.75 ±0.31 
3510.66 ±0.07 

3556.20 ±0.09 

3456.2 

3456.2 

3456.2 

Ip 

Ip 

Ip 

1++ 

2P[To,ri] 

3871.7 

3871.69 ±0.17 

3910.9 

2p 

0++ 

lP[Ti,r2] 

3897.9 

3918.4 ± 1.9 

3910.9 

2p 

2++ 

2P[To,ri] 

3903.0 

3927.2 ±2.6 

3910.9 

2p 

1++ 

lP[Ti,r2] 

4017.3 




0++ 

^P[T3,T4] 

4140.2 

X(4140) 



2++ 

lP[Ti,r2] 

4140.2 




0++ 

^P[T4,n] 

4325.1 

X (4350) 

4294.6 

3p 


Table 5: Calculated charmonium masses from V{r) : Mgspm up to 4350 MeV. The 
(lp[T 2 ,T 3 ]) row has been omitted since there is no CSPM bound state in that energy- 
region. For we do not list any state above 4080 MeV due to the current incomplete 
knowledge about thresholds above this energy. The same for 2^^^ states above 4224 
MeV. Masses for experimental resonances, Mpoc, have been taken from [2] (when a 
resonance appears in the Particle Listing section of [2] but not in the Summary Table we 
write the name of the resonance that contains the nominal mass between parenthesis). 
For p waves we quote separately the npQ, npi and np 2 states. Masses from the Cornell 
potential, Mcor, are also shown for comparison. 
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JPC 


GSPM 

States 

hTi-i,Ti\ 


Mgspm 

MeV 


MpDG 

MeV 


Mcor 

MeV 


Cornell 

States 

k 


l'S[To,Ti] 

3046.0 

3096.916 ±0.011 

3046.0 

Is 

2S[To,Ti] 

3632.2 

3686.09 ±0.04 

3632.2 

2s 


3743.5 

3773.15 ±0.33 

3743.5 

Id 

3S[To,Ti] 

4063.2 

4039 ± 1 

4065.8 

3s 

2'^[To,Ti] 

4139.3 

4191 ±5 

4142.8 

2d 


Table 6: Calculated 1 charmonium masses from V{r) : Mgspm- Masses for experi¬ 
mental resonances, MppQ, have been taken from [2]. Masses from the Cornell potential, 
Mcor, are also shown for comparison. 

(Z)*(2007)°Zi)*°(2007) and D*^D*~) the central potential used does not make any dif¬ 
ference for J = 0 and J = 2. The other new resonance that we shall call (7(4017) 
would be assigned to the {lpyTi,T 2 \) CSPM state at 4017.3 MeV. As shown above 
this resonance is generated altogether with X(3872) as an effect of the introduction of 
the D^D* (2007) threshold. Hence the existence of (7(4017) seems to be unavoidable if 
the mechanism proposed for the generation of X(3872) is the correct one. 

For the sake of completeness let us mention that for energies quite below the hrst 
thresholds the calculated spectrum is of Cornell type giving rise to degenerate = 
(0,1, 2)^"'" states. This degeneracy is broken for energies reaching the hrst thresholds 
(and beyond) due to the different values of the threshold masses in each case. 

We might then conclude that an assignment of CSPM states to the existing well 
established or possible candidates to 0’''(J’''+) resonances is feasible. With respect to 
the observed differences between the calculated CSPM masses and data we shall assume 
that the experimental values can be reached from the CSPM ones through perturbative 
corrections to the hamiltonian. The experimental conhrmation of the candidates and the 
discovery of the new predicted resonances could give dehnite support to this conclusion. 

3.3.2 0(1—) GSPM States 

The spectrum for 0“(1 ) cc states from V{r) up to 4200 MeV (quite below the hrst 

threshold located at 4287 MeV, see Table E]) is shown in Table El The spectrum from 
the Cornell potential Vcorif) given by (jlj) with the same values of the parameters a, y 
and rric is also listed for comparison. 

An almost pure Cornell like spectrum (very little threshold effects) is obtained in this 
energy region as can be checked by comparing the calculated CSPM masses with the 
Cornell ones. 
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GSPM 

States 

Mgspm 

MeV 

Mgot 

MeV 

Cornell 

States 

3'S[ro,ri] 

4063.2 

4065.8 

3s 

2^[To,fi] 

4139.3 

4142.8 

2d 


4337.3 




4454.1 

4436.5 

4s 

2'5[fi,oo[ 

4483.5 

4496.1 

3d 


Table 7: Calculated 1 charmonium masses from 4.0 to 5.0 GeV when only one non 
degenerate threshold Ti = D^Di (2420) is considered: Mgspm- Masses from the Cornell 
potential, Mcor, are also shown for comparison. 


As explained before the GSPM can not be reliably applied to calculate the masses 
of higher spectral states in this case. Nonetheless a qualitative analysis of the possible 
mixing configuration content in some of the well established higher spectral resonances 
can be carried out. Let us centre for instance in X(4260) lying close below the first 
(degenerate) threshold. Let us examine whether this resonance could be obtained or not 
if only the threshold Tn = D^Di^{2420) at 4287 MeV were present. Then the resulting 
GSPM spectrum from 4.0 to 4.5 GeV would be as listed in Table [7] where the Cornell 
spectrum is also given for comparison. 

A glance at the table shows that the presence of the threshold would generate a new 
spectral state at 4337 MeV as compared to the Cornell case. But there would not be any 
chance to obtain a resonance close below threshold, as the X(4260). The explanation 
for this has to do with the fact that there is not any Cornell state close below or above 
threshold from which such resonance could be formed by the effect of the threshold. 

The situation could change by considering the additional effect of the other degener¬ 
ate threshold T 12 = D^Di (2430) on the new spectral state. As the mass of this state, 
4337 MeV, is close above T 12 , it could be shifted down to a value below the threshold 
as experimentally observed (notice that some additional attraction could also be pro¬ 
vided by the Zi)*(2007)°Zi)o°(2400) threshold due to its large width). In this regard a 
refined version of the GSPM, incorporating a lesser abrupt change in the potential when 
approaching the threshold, could allow for a consistent treatment of the degenerate as 
well as the non degenerate threshold effects through the different paths followed by the 
potential to reach the threshold energy. 

Therefore we might tentatively conclude (without any quantitative proof) that the 
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existence of X (4260) could be related to the presence of degenerate overlapping thresh¬ 
olds. Otherwise said X (4260) could be the result of the mixing of the quenched cc with 
D°Di°(2420) and D^Di {243)Q) conhgurations. 


3.3.3 Electromagnetic Widths 

Electromagnetic decay rates of charmonium are sensitive to details of the wave functions 
involved. Therefore their study might serve to test a quark model and to discriminate it 
against others. One should realize though that when ratios of decay rates are considered, 
similar results may be obtained from different models. Indeed some of these ratios, 
involving transitions from initial to hnal charmonium states, can be explained from heavy 
quark symmetry considerations without reference to any particular dynamic model [8]. 

Let us note that the GSPM assigns a differentiated state to each of the existing 
non conventional experimental candidates to be a 0’''( J’*"’*') resonance. It also allows for 
an unambiguous assignment of states to conventional 0+(J’''+) and 0“(1 ) resonances 

below their Erst thresholds. Therefore it can be consistently used for the analysis of 
transitions involving these states. 


We will focus on the calculation of electric dipole (El) and two photon decay widths 
for which a comparative analysis to data can be carried out. Thus, for El decays we 
shall centre on transitions between spin triplet P— wave and S— wave states for which 
the non relativistic El partial widths read [9] 


hiji (z -)■ / -f 7) 


4ae>f^ {2Jf + 1) 
27 


l(/kl*)l^ 


(9) 


where i and / denote the initial (hnal) charmonium state, a stands for the hne structure 
constant, Cc = | is the charm quark electric charge, Wif is the photon energy 

«>./ = ^ - Mj) ( 10 ) 

and (/ |r| i) is the dipole matrix element 

00 

{f\r\i) = J Rf {r) r^Ri (r) dr (11) 

0 


with Rf (r) and Ri (r) standing for the radial wave functions of the hnal and initial state 
respectively. 

From ([9]) we can easily establish the ratios: 
a) 


hgi {i ^ /i + 7 ) ^ '^ffi K/i 

TEi{i^f2 + l) wff,\{f2\r\i)f 
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for the case in which the same initial state decays into two final (/i and / 2 ) states with 
the same value of and 


r^i (^1 ^ / + 7 ) _ ‘^iif \{f M 
rm(^2^/ + 7) wif\{f\r\i2)\^ 


for the case in which two initial states [ii and 12 ) decay into the same hnal state. 


As for two photon transitions we shall consider the decays from and ^P 2 states. 
In the nonrelativistic limit the decay widths can be expressed as 


r (* i^Po) ^ 77 ) = ^ 


r (z (3p2) ^ 77) = 


m- 

36a^e^ 

5mt 


R. 


R. 


i^Po) 


(3P2) 


( 0 ) 


( 0 ) 


(14) 

(15) 


where R[ (0) stands for the derivative of the radial wave function at the origin. 

First order QCD radiative corrections to (1141) and ([IS]), in the form of multiplying 
factors, have been calculated. For the effective value of as = 0.38 we are using they are 
signihcant. This poses the need to calculate them to higher order. Instead we shall keep 
the zeroth order expressions to get a hrst approach to data and we shall use for practical 
purposes the ratios 


r {t2 i^Po) ^ 77 ) 

F(^i (3 Po)^77) 

r {i 2 {^P 2 ) 77 ) 

r(*i (3P2)^77) 

where the multiplying factors cancel out. 



R' 


n(3P2) 


( 0 ) 


(16) 


(17) 


\'oj(lP) 

The Xcj i^P) resonances are identified with the (lP[ro,ri]) GSPM states which 
are practically identical to the {Ip) Cornell states. As the J/zp description is also 
the same with both models they give the same results for the XcJ (1-P) —)■ decay 

widths if the same values for the photon energies are chosen. 

The calculated GSPM widths are shown in Table [H] where the dipole matrix elements 
are also tabulated. For the photon energies the experimental values have been used. 
This can be justihed under our former assumption that the experimental masses can 
be reached from the GSPM ones through hrst order perturbative corrections to the 
hamiltonian (let us remind that no modihcation of the wave functions is then generated). 
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/ 

i 

{'^if) Exp 

MeV 

fm^ 

(Pei (^ / -1 - x))gspm 

MeV 

(Piji {i ^/ + 7)) 

MeV 

J/^ 

XcO (IP) 

303.04 

0.198 

0.20 

0.13 ±0.02 


Xcl (IP) 

389.36 

0.198 

0.43 

0.29 ±0.02 


Xc2 (IP) 

429.63 

0.198 

0.58 

0.37 ±0.04 


Table 8: Calculated El dipole matrix elements (fourth column) and decay widths (hfth 
column) for Xcj (1-P) yJ/'i/'- Photon energies (third column) from ffTOj) with experi¬ 

mental masses |2]. Widths data (sixth column) from [2]. 


As can be checked the values obtained are 30% off the experimental intervals. This 
can be considered a reasonable hrst approach to data and a starting point to include 
additional corrections (see [13] and references therein). Although we do not proceed here 
along this line it is worth to point out that the central values of the experimental ratios 


(Pei iXcoilP) ^ 1J/'>P))Exp 

and 

(Fe, (x=. (IP) ^ „ 

are in good agreement with the experimental photon energy ratios 


(18) 


(19) 


w. 


Xc2(\P)J jiii 


Exp 


W. 


XcO(iP)J/i> 


Exp 


2.9 


(w^ ) ~ ' 

V Exp 

Taking into account flT^ this suggests the additional corrections should not introduce 
any signihcant difference among the Xco (1-P), Xci (1-P) and Xc 2 (IP*) wave functions. 


Regarding two photon decay widths, the degeneracy of the J++ (lp['ro,ri]) GSPM 
states gives rise from fllTll and flT5l) to the ratio 


(r (Xco (1-P) h))gspm _ 

(P {Xc2 {IP) ^ 11))gspm (f ) 
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Figure 3: Radial wave functions R(r) (in units fm 2 ) for the 2++ GSPM state 

(thick line) and the {2p) Cornell state (thin line). 


within the experimental interval [2] 

(F iXcoilP) ^77 ))exp ^ 2.3 ± 0.4 KeV 
(r (Xc 2 (1^) ^ 11 ))Exp 0-53 ± 0.07 KeV 


Xc2 (2P) 

The Xc 2 (2P) resonance at 3927 MeV is assigned to the 2++ (2p[To,Ti]) GSPM state 
that differs little from the 2'^^ {2p) Cornell one, as shown in Fig. [3l 

Only an experimental lower bound for the two photon decay width is known |2] 

(F (Xc 2 (2p) ^ 11))Exp > 0-21 ± 0.04 eV 

from 

(F (Xc2 (2p) ^ 77) B {Xc 2 (2p) ^ = 0.21 ± 0.04 eV 

where B stands for branching fraction. 

From the calculated GSPM wave functions we get from flTTll the ratio 

(r (xc 2 (2p) ^ 77)) GSPM 
(r (Xc2 (Ip) ^ 77)) GSPM 


^2++{2p^r„T,{) 


GSPM 


= 1.34 


^( 0 ) 

2++(lp[To,Ti]) ^ 


GSPM 
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By assuming that this value is a reasonable approach to the experimental ratio we might 
expect the approximated values 

r (Xc2 (2p) ^ 77 ) - 1-34 (r {Xc2 (Ip) ^ 'll)) Exp = ± 0-09 keV 

B {xc2 (2p) ^ DD) ~ 0.30 ± 0.10 


X(3872) 

As shown before the I’*""'' (2p[To,Ti]) GSPM state is identihed with the X(3872) whose 
mass has been used to hne tune the charm quark mass. 

Concerning electromagnetic decays the ratio 

r(X(3872)^7^(2g)) 
r (X(3872) ^ xJ/^) 
has been recently measured |T5] to be 

Aexp = 2.46 ± 0.64 ± 0.29 

compatible with the previous value 3.4 ± 1.4 |Tn] and the upper bound < 2.1 mi. 

From flT^ the GSPM gives for this ratio the value 

Agspm = 2.01 

calculated from the dipole matrix elements 

KV’ (2s) |r| X(3872))ggp^|^ = 0.2856 W 
\{J/^P\r\X{m2))^SPM? = 0-0025 W 

and the experimental values of the photon energies 

('?i'x( 3872 )p( 2 s))p^p = 181.25 MeV 
(wx( 3872 )j/p)p^p = 697.19 MeV 

Therefore a full compatibility with existing data comes out. We should point out though 
that quite the same result would be obtained for the dipole matrix elements by using 
the !+■'■ (2p) Cornell state wave function instead of the l’''+ {2px:'o,Ti]) GSPM one. As 
the main difference between these two wave functions is the long tail of 2p]^To,Ti] state 
as compared to that of 2p (see Fig. |2]) we may conclude that these radiative decays are 
not sensitive to the long distance nature of X(3872). This can be understood by the 
negligible long distance overlap of the 1+’'' (2p['7|,-rp) state with J/ip and ip (2s). (The 
calculated root mean square radii for J/ip and ip (2s) are respectively 0.5 fm and 0.9 fm.) 
It should be mentioned that the same conclusion has been also inferred by other authors 
using molecular descriptions for X(3872) [T8l 119] . 
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Figure 4: Radial wave functions R(r) (in units fm 2 ) for the GSPM state 

(thick line) and the {2p) Cornell state (thin line). 


X(3915) 


The Review of Particle Properties [2] has identihed the X(3915) with a conventional 
Xco (2p), this is with a O’*"’*' (2p) Cornell like state. This identihcation has been criticized 
by some authors [201 EH- A major criticism is the lack of evidence of X(3915) —?■ DD 
decays. From our estimation above for the branching fraction B (xc 2 (2p) —)■ DD^ ~ 0.30 
a similar result could be expected for B (xco (2p) —!■ DD) since the 0+’'' (2p) and 2++ (2p) 
Cornell states are degenerate and the measured values of the masses and total widths of 
Xc 2 (2p) and X(3915) are quite similar. 


m 


Regarding electromagnetic processes a lower bound for the two photon decay width 

(r(X( 3915 )^ 77 ))B.p> 54 ± 9 eV 


is known from 


(F (X(3915) ^ 77 ) B (X(3915) ^ J/^oj))exp > 54 ± 9 eV 


The CSPM assigns the X(3915) to the O’*""'' (lp[Ti,r 2 ]) state that differs greatly from 
the O’*"’*' (2p) Cornell one as shown in Fig. 01 

From the calculated CSPM wave functions we get from flTB]) the ratio 


(F (X(3915) ^ 77)) GSPM 

{r(xa(iP)^n)) GSPM 


^( 0 ) 

0++(lp[Ti,T2l) ^ ^ 

4(0) 

0++(lp[To,Ti]) ^ 


0.02 
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Assuming again that this value is a reasonable approach to the experimental ratio we 
might expect 


r (X(3915) -)■ 77 ) ~ 0.02 (r {Xco (Ip) ^ 11))Exp = 44 ± 7 eV 
By combining this result with the experimental lower bound given above we would get 

r (X(3915) ^ 77 ) ~ 48 ± 3 eV 
i3(X(3915) ^ J/iJu) > 0.88 

It should be emphasized that the identihcation of X(3915) with O’*""'' (2p) would give 
a completely different ratio 

(r (X(3915) ^ 77 )) Cornell 

(r (y,o(lP)^77)) 

Cornell 

and consequently completely different values for the two photon decay width (2640 eV) 
and the branching fraction to J/tjjuo (> 0.02). 

Certainly the big GSPM branching fraction for X(3915) —)■ J/ipoj, OZI suppressed in 
the Cornell model, should be somehow justihed. A quantitative justihcation, if possible, 
would imply the development of a strong decay theory above threshold within the GSPM 
framework which is outside the scope of this work. Hence we shall limit here to a merely 
speculative qualitative comment. Let us imagine for instance that due to the threshold 
modulation the GSPM interaction favored, at the energy of X(3915) and through light 
quark pair creation out of the vacuum, the formation of color octets made of heavy 
quark-light antiquark and viceversa. Then it would be possible to have a dominant decay 
through reordering of the quarks in such a state like On the contrary the Cornell 

states are known to favor the formation of color singlets giving rise to the dominant fall 
apart decay mode DD. 

Therefore the GSPM and Cornell descriptions represent incompatible scenarios for 
the understanding of X(3915). According to our analysis more detailed data could deh- 
nitely clarify the situation about the true nature (non conventional or conventional) of 
this resonance. 




0++(2p) 


( 0 ) 




0++(lp) 


( 0 ) 


= 1.2 


4 Summary 

A nonrelativistic quark model called Generalized Screened Potential Model, or abbreviate 
GSPM, previously used to calculate the bottomonium spectrum has been applied to 
charmonium. 
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The model, whose interaction potential has a Cornell form but modulated by meson- 
meson thresholds, has been used to calculate 0’''( charmonium masses up to 4.4 GeV, 
a limit imposed to the application of the model by the incomplete current knowledge of 
open charm meson-meson thresholds. As it turned out to be the case in bottomonium 
a richer spectrum (bigger number of bound states) than the one resulting from the non- 
screened Cornell potential is predicted. However, differing from bottomonium where 
the lack of data prevented the verihcation or refutationof such a prediction, there exist 
in charmonium well established as well as candidates to non conventional resonances 
in the energy interval analyzed. As a matter of fact the well established X(3872) is 
nicely described as a GSPM state that can be interpreted as being generated from the 
D^D* (2007) threshold and the (2p) Cornell state. Regarding the experimental 
candidates X (4140) O’*" and X (4350) O’*" (7^''') a good spectral correspondence 
with GSPM states is observed. Furthermore two new 0’''( J’*"’*') resonances are predicted, 
a 0’'“(2’''’'') or 0’'“(0’'“’'') one with mass around 4140 MeV and a 0’''(1’'“’'') one with mass 
around 4017 MeV. (Notice though that all the calculated GSPM masses, except for 
A(3872) which is used to hne tune the quark mass, are below the experimental ones 
what suggests that the masses of these new resonances could also be underestimated.) 
The generation of these new resonances in the GSPM is related to the presence of the 
states assigned to X(4140) and X(3872). Therefore their discovery would constitute a 
dehnite check of the GSPM as a model for the spectral description. 

The GSPM has also been employed to evaluate the 0“(1 ) spectrum up to 4.2 GeV, 

a limit imposed to the applicability of the model by the presence of degenerate and 
overlapping thresholds. The resulting spectrum in this energy region is very much of 
Gornell type. In order to go further in energy the model should be rehned. With respect 
to this a simplihed qualitative analysis of the possible generation of X (4260) seems to 
point out that the joint effect from overlapping meson-meson conhgurations should be 
an essential ingredient to be incorporated. 

A study of electromagnetic decays of 0'''( J'*""'') resonances, specihcally El and two 
photon decays for which there are data available, has also been carried out. The calcu¬ 
lated GSPM widths are fully compatible with existing data. However, more detailed data 
are needed to perform a stringent check of the GSPM. In this regard a thorough exper¬ 
imental analysis of X(3915) is of particular interest given the very different description 
coming out from the GSPM and the Gornell models. 

This work has been supported by Ministerio de Economia y Gompetitividad of Spain 
(MINEGO) grant FPA2013-47443-G2-1-P, and by PrometeoH/2014/066 from Generali- 
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